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ABSTRACT 

Many  separated  targets,  of  only  a  few  differing  values,  are  subject  to  a  simultaneous  attack. 
The  area  defenses  considered  have  (a)  impact  point  prediction  (IPP)  and  full  coordination,  or  (b)  no 
IPP  and  full  coordination,  or  (c)  no  IPP  and  partial  coordination.  For  a  given  attack,  the  defense 
wishes  to  allocate  its  interceptors  to  maximize  the  expected  total  survival  value  of  the  targets.  For  a 
given  attack  size,  and  with  knowledge  of  the  defense’s  capabilities,  the  offense  seeks  a  strategy  to 
minimize  expected  total  survival  value  against  best  defense.  We  present  algorithms  to  determine 
optimal  a' tack  and  defense  strategies  and  the  optimal  value  of  the  min-max  problem,  and  we  show 
how  to  take  computational  advantage  of  the  relatively  few  unique  target  values.  Illustrative 
computational  results  are  provided.  ■  .  ■  ■  L  ■  ■  ■ 


1.  INTRODUCTION 


Consider  an  area  to  be  defended  that  consists  of  T  separated  point  targets,  and  associate  with 

target  i  a  value  v,,  i  =  1,  T.  Without  loss  of  generality,  we  assume  that  the  v,  satisfy  vt  >  v.,  > 

■  >  v-p  >0.  We  denote  the  number  of  attackers  directed  against  target  i  by  a,  and  the  attack 

T 

against  the  T  targets  by  a=  [a^  aT],  with  ^  a,  =  A  and  each  a,  e  K,  the  set  of  nonnegative 

i=l 

integers.  The  defense  has  D  interceptors,  each  of  which  may  be  assigned  to  intercept  any  attacker  in  a 
given  subset,  possibly  the  entire  set,  of  attackers.  The  objective  of  the  defense  is  to  allocate  its 
interceptors  against  attackers  so  as  to  iiiaxm.izc  expected  total  survival  value.  In  his  first  move,  the 
attacker  seeks  an  attack  allocation  which  will  minimize  expected  total  survival  value  against  best 
defense. 


The  specific  assumptions  we  make  are  as  follows: 

(1)  The  attack  is  simultaneous,  i.e.,  the  defense  sees  the  overall  size  of  the  attack,  A. 

(2)  The  defense  has  the  “last  move.’'  Only  after  the  attack  is  observed  does  the  defense 
need  to  allocate  its  D  interceptors. 

(3)  All  A  attackers  are  identical.  An  attacker  will  penetrate  to  its  intended  target  if 
not  successfully  intercepted.  Such  a  penetrator  will  destroy  the  target  with  probability  a,  a  £  (0,  1]. 

If  multiple  attackers  penetrate  to  a  target,  they  act  independently. 

(4)  All  D  defenders  (interceptors)  are  identical.  An  interceptor  will  destroy  the 
attacker  at  which  it  is  directed  with  probability  6,  6  £  (0,  1].  If  several  defenders  are  aimed  at  a  single 
attacker,  they  act  independently. 


□  □ 


(8)  If  the  defense  has  impact  point  prediction  (IPP),  it  correctly  identifies  the  intended 
target  of  each  attacker. 

(9)  If  the  defense  is  (fully)  coordinated,  its  resources  are  allocated  by  a  single  controller. 
If  different  subsets  of  the  interceptors  are  allocated  by  independent  controllers,  we  say  the  defense  is 
partially  coordinated. 

(10)  The  defense  has  area  defenders  only;  none  of  the  defensive  interceptors  are  point 
defenders  assigned  to  protect  specific  targets. 

A  brief  outline  of  the  paper  is  as  follows.  In  the  remainder  of  this  introductory  section  we 
briefly  survey  the  literature  dealing  with  problems  of  attack  and  defense  of  a  set  of  separated  point 
targets.  In  section  2,  we  introduce  our  problem  functions  and  formulate  the  defender’s  problem  and 
the  attacker’s  problem  associated  with  each  of  the  three  defense  capabilities  considered:  (1)  the  defense 
has  IPP  and  is  coordinated;  (2)  the  defense  has  no  IPP  but  is  coordinated;  and  (3)  the  defense  has  no 
IPP  and  is  only  partially  coordinated.  For  each  of  the  three  cases,  we  provide  (or  reference) 
appropriate  algorithms  to  solve  both  the  defender’s  problem  and  the  attacker’s  problem  and  provide 
the  optimal  min-max  value.  In  section  3  we  show  how  the  efficiency  of  each  of  the  three  algorithms, 
for  the  respective  versions  of  the  attacker’s  problem,  can  be  greatly  increased  by  taking  advantage  of 
the  relatively  few  distinct  values  among  the  T  target  values.  Section  4  provides  illustrative 
computational  results,  and  section  5  contains  concluding  remarks. 

A  good  deal  of  missile  allocation  work  deals  with  preallocated  preferential  defense  of  a  set  of 
identically  valued  targets,  using  a  game  theoretic  approach  by  both  attacker  and  defender.  Early  work 
on  this  problem  was  by  Matheson  (1966,  1967)  and  recently  Bracken,  Brooks  and  Falk  (1987)  and 
Bracken,  Falk  and  Tai  (1987)  have  examined  robustness  issues  connected  with  the  assumption  that  the 
defender  knows  the  attack  size  A.  They  provide  references  to  the  intervening  work  in  this  area. 

In  their  seminal  work,  Eckler  and  Burr  (1972)  examined  a  number  of  missile  allocation 
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problems  under  a  variety  of  assumptions  and  thoroughly  surveyed  the  literature. 


For  targets  of  different  value,  Karr  (1981)  analyzed  several  problems  under  the  assumption  of 
sequential  attack  (the  defense  sees  the  attackers  one  at  a  time  and  must  make  an  engagement  decision 
at  that  time  without  knowledge  of  the  overall  attack  size)  and  developed  continuous  (i.e.,  noninteger) 
solutions  on  the  basis  of  Prim-Read  deployments. 

There  is  a  body  of  work  in  which  the  defender’s  objective  is  to  minimize  the  cost  of  the  defense 
while  insuring  prescribed  levels  of  survival  value  as  a  function  of  attack  size.  Under  the  assumption  of 
a  linear  bounding  function,  Burr,  Falk  and  Karr  (1985)  generated  integer  solutions  for  a  defense  against 
a  sequential  attack.  Falk  (1985b)  extended  this  concept  to  a  larger  class  of  bounding  functions,  and 
McGarvey  (1987)  added  an  overview  with  mixed  strategies.  Karr  (1982)  and  Soland  (1987a)  developed 
solutions  for  defenses  against  simultaneous  attacks. 

Finn  (1986)  examined  the  value  of  the  additional  knowledge  provided  by  IPP  for  the  case  of 
perfect  weapons  (i.e.,  a  =  6  =  1)  with  an  area  defense  and  terminal  defenses  at  the  targets. 

2.  PROBLEM  FORMULATIONS  AND  SOLUTION  ALGORITHMS 

2.1  Coordinated  Defense  with  IPP 

2.1.1  Problem  Formulation 

If  d  •  interceptors  are  allocated  for  the  defense  of  target  i  under  an  attack  of  size  a-,  then  an 
even-as-possible  distribution  of  these  d(-  defenders  against  the  a-  attackers  is  optimal.  [That  such  a 
quasi-uniform  defense  of  target  i  is  optimal  is  shown  by  Soland  (1987b)].  Specifically,  for  a-  >  0,  the 
defense  allocates 

fd  -  /a;l  interceptors  to  each  of  =  [df  —  [dj/a^jaj  attackers 

and 

k,  =  [dj/aj  interceptors  to  each  of  a(  —  r;  attackers. 

[Notation:  We  use  [xj  for  the  integer  part  of  x  and  [x]  =  min  {n  >  x  |  n  integer}.] 
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If  we  define  the  probability  that  a  target  survives  an  attacker  facing  m  defenders  as 


iTrr.  —  1  —  a(I—  6)m  for  m  =  0,  1,  ...  , 

then,  with  k(  and  rt  as  defined  above,  we  may  write  the  survival  probability  at  target  i  as 

f  a,  =  1,2,  ...,  A 

d<)  =  { 

l  1,  a,  =  0. 


Let  d  =  [dj,  ...»  d-f].  We  now  define  the  defender’s  problem  for  a  given  attack  a  as 

T 

?x(a)  Maximize  $>(  a,  d  )  =  £  v,0(a,,  d,) 

i=i 

subject  to 

E  d,  <  D 


di  £  X,  *=1,  ...  ,  T. 


We  denote  the  optimal  value  of  problem  ^i(a)  by  v(9j(a)].  [Note:  For  perfect  weapons  (a  =  6  =  1), 
the  defender’s  problem  is  a  0-1  knapsack  problem.  Otherwise,  it  is  a  nonlinear  integer  programming 
problem.]  We  may  then  formulate  the  attacker’s  problem  as 

Minimize  v[  ^j(a)] 
subject  to 

Ea<  <  A 

i=i 


a<  £  K,  t  =  1,  ,..,T. 


Inasmuch  as  our  solution  algorithm  for  problem  ?!  (itself  a  nonlinear  integer  programming  problem)  is 
of  the  branch-and-bound  type,  we  would  greatly  prefer  to  limit  the  number  of  attacks  that  we  need  to 
examine.  The  following  Nonincreasing  Attack  Theorem  (NAT)  enables  us  to  do  so.  The  NAT  for  the 
case  of  IPP  and  perfect  weapons  was  first  proved  by  Falk  (1985a),  whereas  a  proof  with  IPP  and 
imperfect  weapons  (i.e.,  a  <  1  and  6  <  1)  was  suggested  by  Falk  (1987)  and  given  by  O’Meara 
(1988).  The  proof  for  the  case  of  no  IPP  and  coordinated  defense  was  given  by  O’Meara  and  Soland 
(1988b). 
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Theorem  (Nonincreasing  Attack  Theorem) 


There  is  an  optimal  attack,  say  a*,  with  a*  >  a*+1,  »  =  1,  ....  T  —  1.  We  call  a*  a 
nonincreasing  attack. 

2.1.2  Algorithm  for  Problem 

The  algorithm  we  use  to  solve  problem  is  a  branch-and-bound  scheme.  During  the  course 
of  the  algorithm,  it  is  often  necessary  to  determine  the  “best  defense”  against  some  given  attack  a,  that 
is,  to  solve  problem  ^i(a)  (or  to  determine  some  bounds  on  its  optimal  value).  Details  of  the 
algorithm  are  provided  by  O’Meara  and  Soland  (1988a,  1988c). 

2.2  Coordinated  Defense  with  no  IPP 
2.2.1  Problem  Formulation 

To  avoid  any  cumbersome  notation,  we  shall  examine  the  cases  of  perfect  weapons  (ar  =  l,6  =  l) 
and  imperfect  w  apons  (a<l,  6<1)  separately.  Results  for  the  mixed  cases  follow  without  difficulty. 

From  the  point  of  view  of  the  attack,  with  perfect  weapons  there  is  a  logical  constraint  of 
allocating  no  more  than  D  +  1  attackers  against  any  target  since  that  number  will  assure  target 
destruction.  The  assumption  of  no  IPP  means  that  the  defense  has  no  real  decision  since  all  the 
attackers  “look  the  same.”  Accordingly,  the  defense  allocates  its  D  defenders  against  D  of  the 
attackers.  Since  the  weapons  are  perfect,  ( A  —  D )  attackers  will  penetrate  to  their  intended  targets. 
For  nontrivial  problems  we  have  D+l  <  A  <  T(D  +  1). 

We  observe  that  target  i  survives  only  if  all  a-  attackers  are  intercepted.  To  determine  this 
probability,  we  consider  the  D  attackers  to  be  intercepted  as  a  sample  drawn  at  random  (without 
replacement)  from  the  population  of  A  attackers.  If  we  think  of  this  as  a  population  dichotomized  into 
attackers  destined  for  target  i  and  those  that  are  not,  a  hypergeometric  distribution  is  appropriate. 
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Accordingly,  we  define  the  survival  probability  at  target  «  against  attack  <r  as  ip(a,)  and  compute  it  as 


^  0,  a<  >  D. 

We  then  compute  the  expected  survival  value  at  target  i  as 

v.V>(a  i) 

and  the  expected  total  survival  value  for  attack  a  as 

#(*)  =  £v^(ai)- 
»=i 

With  imperfect  weapons,  any  number  of  attackers  may  be  directed  against  any  target,  and  any 
number  of  interceptors  may  be  allocated  against  any  attacker.  Of  course,  the  overall  attack  and 
defense  size  constraints  are  still  active. 

As  above,  the  lack  of  impact  point  prediction  means  that  the  attackers  appear  identical  to  the 
defense.  Intuition  suggests  that  the  best  defense  is  to  spread  the  interceptors  as  evenly  as  possible 
among  the  attackers.  This  has  been  shown  by  Soland  (1987b)  to  be  an  optimal  defense  in  that  it 
stochastically  minimizes  the  number  of  attackers  that  penetrate  the  defense.  Specifically,  the  defense 
allocates 

[D/ A  ]  interceptors  to  each  of  R  =  [D  —  [D/A  j  A  ]  attackers 

and 

[D/A  J  interceptors  to  each  of  A  —  R  attackers. 

With  system  reliabilities  a,  6  z  (0,  1),  the  parameter  range  for  nontrivial  problems  is  A  >  0  and 
D  >  0. 

In  this  case  the  uncertainty  about  survival  value  can  be  attributed  to  both  system 
(un)reliability  and  the  lack  of  IPP.  As  might  be  expected,  the  survival  probability  function  has  a  more 
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involved  structure.  For  target  i  not  under  attack,  we  have  V'(»,)  =  1-  For  a, 
arguments  (see  O’Meara  and  Soland  (1988b)]  yield  the  following: 

(i)  if  I)  <  A, 


a,) 


/  —  n  n  =  n '  *  n  ' 


where  rj  =  max  {a,  —  D,  0  }  and  ti  =  min  {/,  A  —  D  }; 


(ii)  if  D  >  A, 


V’(a>)  =E(1-f>)1 

f=n 


■t  £U)  Lu -*) 

m  =  0  n  =  n 


— .m  r 


,at -n-(  t-,-11 


fix* 

(a-r) 


wuert  n  —  max  {a;  —  R,  0  }  and  ?i  =  min  {a,  —  (/  —  m),  A  —  R }. 


We  formally  define  the  attacker’s  problem  as 


^2 


T 

Minimize  'l'(a)  =  ^v^aj 
i- 1 

subject  to 

2>.  <  A 

i=  l 


a,  e  Jf,  T. 


0,  conditioning 

). 
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While  it  is  fao.y  straightforward  to  show  that,  for  perfect  weapons,  problem  (5P2  is  a  convex 
programming  p.oblem  in  the  integer  sense  that  the  function  4'  has  nondecreasing  increments,  the 
complexity  of  4>  for  imperfect  weapons  has  not  yet  permitted  a  similar  theoretical  assertion.  However, 
all  computational  experience  thus  far  has  not  uncovered  a  nonconvex  case.  The  algorithm  we  use  (see 
below)  is  of  the  “reduction”  type  for  convex  problems.  A  dynamic  programming  algorithm  is  also 
available  for  any  nonconvex  instances,  if  necessary.  [See  O'Meara  and  Soland  (1988b)]. 

2.2.2  Algorithm  for  Problem  ^P2  (V’  Convex) 

Lafon  and  Lahrichi  (1988)  describe  the  following  algorithm  for  a  class  of  problems  cont aiuing 
problem  when  V  has  nondecreasing  increments,  that  is,  when 

V(a,  +  1)  -  v(a,)  >  t/(at)  -  V’(a.-l),  a,  =  1,  2 . A. 

Algorithm  [Lafon  and  Lahrichi  (1988)]: 

Step  0  (Initialization): 

a,  =  0,  i  =  l . T:  A,  =  A; 

Jt=l;  go  to  step  k. 

Step  k  ( Allocat ion): 

If  A t  =  0,  stop  (a  is  optimal  for  problem  ^P2); 
else  xt  =  r Afc/T]; 

i*  =  minji  |  v,[t'(a,  +  xt)  -  tfa,)]  =  min  y  {vy[t’(a;  +x* )  -  f(ay)]}.  i=l . t|: 

a,.  :=  a,-  +  xt; 

^t+i  :=  'V+i  _  xkl 
k  :=  k+  1;  go  to  step  k. 
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Ohs.  •rve  that  the  allocation  a  =  [a, ,  ....  a-j-]  at  the  start  of  any  step  Ic  >  1  solves 


T 

Minimize  ]T  v,  i/’(a,) 
subject  to 

I>.  —  A  —  Ak 
1  =  1 

a,  £  Jtr,  t=l . T. 


Partially  ( 'oordinatcd  Defense  with  no  1PP 
-.•1. 1  Problem  Formulation 

For  the  case  of  partial  coordination  and  no  IPP,  suppose  the  D  interceptors  are  divided  into  N 
(b'lltlr)  group s  of  n  interceptors  each,  so  Nil  =  F).  Suppose  each  group  of  n  interceptors  is  fully 
coordinated  in  the  sense  that  it  is  centrally  controlled.  We  assume,  however,  that  different  groups  act 
independently  and  do  not  communicate  with  each  other. 

Suppose  each  of  the  N  groups  is  physically  able  to  engage  each  of  nt  different  attackers  with 
each  of  its  ri  interceptors,  where  1  <  m  <  A. 

Assume,  finally,  that  each  attacker  may  be  defended  against  by  r  different  groups,  and  thus  by 
a  total  of  nr  interceptors.  Hence,  there  are  Anr  possible  pairings  of  attackers  and  interceptors.  This 
number  is  also  equal  to  mD.  since  each  interceptor  can  be  allocated  to  any  one  of  m  different  attackers. 
Hence  Anr  =  mD,  so  n/rn  =  D/rA.  For  example,  A  =  D  =  1000,  n=10  and  m=20  lead  to  r  =  2. 

Now  consider  an  arbitrary  one  of  the  a,  attackers  directed  at  target  a,  and  consider  one  of  the  r 
groups  of  interceptors  able  to  fire  at  this  attacker.  Since  the  defense  lacks  IPP,  the  controller  of  this 
group  can  do  no  better  than  to  spread  iiis  n  interceptors  as  uniformly  as  possible  among  the  m  attack¬ 
ers  he  can  fire  at;  this  stochastically  minimizes  [see  Soland  (1987b)]  the  number  of  attackers  (out  of  the 
m)  that,  in  the  eyes  of  this  controller,  penetrate  his  defense.  Thus  the  controller  selects  n  —  m[n/mj 
of  the  m  attackers  to  receive  ln/m]  =  [a;]  interceptors  each,  where  t]  =  n/m  =  D/rA,  and 
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m  —  n  -f  m[n/mj  of  the  m  attackers  to  receive  interceptors  each.  Our  particular  attacker  there¬ 
fore  survives  the  defense  of  this  group  with  probability 


p  =  (1  —  6)^\n  —  m[n/mj)/m  +  (1  —  —  n  +  rr>Ln/mJ)/m 

=  (1  -  6){rl\ri  -  L»?J)  +  (1  -  6)[r?J(i  -  n  +  L*?J), 


and  it  survives  the  defenses  of  all  r  groups  that  can  fire  upon  it  with  probability  pr  because  the  groups 
act  independently. 


The  probability  our  particular  attacker  fails  to  destroy  target  i  is  then 

7=1-  apr, 

where,  it  turns  out,  7  <  1  if  6  <  1  and/or  D/rA  <  1. 

It  is  certainly  not  true  that  the  a,-  events  {attacker  j  fails  to  destroy  target  i } .  j  =  1.  ....  a,,  are 
independent.  However,  if  the  a.i  attackers  are  cross-targeted  (i.e.,  launched  from  d’ffercnt  places),  and 
if  a,/A  is  small,  then,  to  a  very  good  approximation,  these  a,  events  may  be  treatec  a  independent. 
This  is  akin  to  approximating  sampling  without  replacement  by  sampling  with  replacement. 


Using  the  approximation  just  described,  we  conclude  that  target  1  survives  with  probability 
7  '.  The  total  expected  survival  value  for  attack  a  is  then  T(a)  =  £  v,7  *,  so  that  the  attacker’s 
problem  against  a  defense  with  partial  coordination  is  formulated  as 


<3>3  Minimize  T(a)  =  v,  1*' 

i=l 

subject  to 

£a.  <  A 

1  =  1 

a,  e  K,  1=  1,  ...,  T. 

For  7  e  (0,  1)  the  function  7  *  for  a<  e  N  is  easily  shown  to  be  strictly  decreasing  with  strictly 
increasing  increments.  This  is  sufficient  to  conclude  that  the  NAT  holds  for  problem  ^P3. 
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Note  that  r,  rather  than  N,  is  the  appropriate  measure  of  the  defense’s  partial  coordination  and 
that  in,  and  hence  r,  depend  on  geometrical  considerations.  Since  73'  is  a  decreasing  function  of  r, 
larger  values  of  r  indicate  less  defensive  coordination;  r=l  indicates  perfect  defensive  coordination  and 
yields  solutions  to  *?3  that  are  almost  identical  to  the  corresponding  solutions  to  <3>2-  The  range  of 
meaningful  values  of  r  is  1  <  r  <  D. 

2.3.2  Algorithm  for  Problem  ^P3 

The  principal  motivation  behind  the  algorithm  for  Problem  1P3  is  that  we  can  take  advantage 
of  the  exponential  property  of  the  survival  probability  function.  The  change  in  expected  survival  value 
at  target  i  when  the  attack  changes  from  a;  to  a{  +  A  is  v,  —  1),  A  =  1,2,....  The  algo¬ 

rithm  is  of  the  “greedy”  type.  That  is.  the  allocation  of  each  successive  attacker  is  made  to  that  target 
offering  the  greatest  immediate  decrease  in  expected  survival  value,  given  the  allocation  to  that  point. 

Algorithm: 

Step  1  (Initial  Phase): 

a,  =  0,  i=l . T; 

determine  T'  =  max  {/  |  £  [  I  +  £n(v,/v,)/Cn7  )J  <  a|; 

,  =  1 T  *•  i  =  l  t  J 

a,  =  [  1  +  fn( vt/v, )/f  117  J,  i  =  1 . T';  A'  =  A  -  £  a,; 

1  =  1 

aT,  =  minfl,  A7};  A1  :  =  A1  —  aT,; 

if  A  —  0,  stop  (a  is  optimal  for  jP3);  else  go  to  step  2. 

Step  2  (Intermediate  Phase): 

a,  =  a,  +  |_a'/T'J,  i=1 . T';  A'  :=  A1-  T'IA'/T'J; 

if  A1  =  0,  stop  (a  is  optimal  for  1P3);  else  go  to  step  3. 
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Step  3  (Final  Phase): 


Let  »,,  ij,  he  a  reordering  of  1,  such  that 

if  1  <  m  <  u  <  Tf,  then  vm7a"*  >  v„7a,>; 
at  -  —  ak  F  1 !  k — 1,  •  ••, 

stop  (a  is  optimal  for  <53). 

3.  GROUPS  OF  EQUAL-VALUED  TARGETS 

If  there  are  relatively  few  distinct  values  among  the  target  values  v,,  v-p,  the  algorithms  of 
section  2  can  be  modified  to  run  much  faster.  The  basis  for  the  increase  in  speed  is  the  following 
Quasiumfonn  Attack  Theorem  (QAT). 

Theorem  (Quasiuniform  Attack  Theorem) 

With  an  appropriate  survival  function  (including  those  of  problems  *3*,,  ^P2  and  ^P3)  and 
equal-valued  targets,  there  is  a  nonincreasing  optimal  attack  a*  with  the  property  0  <  a*  —  a^-  <  1. 

Proofs  of  the  QAT  for  problems  9,  and  are  given  by  O’Meara  (1988).  The  proof  for 
problem  ^P3  follows  directly  from  the  convexity  of  the  survival  probability  function  . 

3.1  Modification  of  the  Algorithm  for  Problem 

To  make  the  necessary  branching  modification  to  the  algorithm  described  by  O'Meara  and  So- 
land  (1988a,  1988c),  we  introduce  the  following  notation.  Order  the  G  different  values  among  v,,...,  v-p 
as  t),  >  •  •  •  >  vg  >  ■  ■  -  >  vG 

with  n$  targets  having  value  vg,  g  =  1,  ...  ,  G.  Furthermore,  define 

tj  -  { *lv,  =  vg,  i  =  1,  ...,  T}, 

Ag  =  a(  and  A9  =  £  *c 
|C  r g  icT g 
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WIh'ii  a  given  node  S  =  {a|  a  ■<  a  ■<  a }  is  presented  for  branching,  we  now  define  the  index 
g*  =  minj<7|0  <  Ag  -  Ag  =  min{A;  -  A;  >  0|>  e  {1,  G}}|  . 


We  then  define  the  limiting  aggregate  values  of  S*~  as 


Air  =  Atf  V  g  =  1 . G  and  A^  = 


and  the  limiting  aggregate  values  of  SR  as 


L(A,  +  A,)/2J,  g  =  g* 
Aj,  g  #  g* 


AR  =  A,  V  g  =  1,  G  and  AR 


L(A,  +  A#)/2J  +  1,  g  =  g* 
As,  3  ^  g‘- 


To  maintain  nonincreasing  limiting  attacks,  we  tighten  according  to  the  following  rules.  In  succession, 
for  g  —  g* . G,  we  tighten 


1 

r 

>1 

»cr 

:=  mini 

1  al 

^9  » 

n« 

Aj-i 

L  n*->  J 

) 

and  for  g  —  g*—  1 .  1,  we  tighten 


D 

A  g  :  —  max 


0-1  _  _ 

Define  N,  =  0  and  Ny  =  n,  for  g  =  2,  ...,  G.  Then  with  tg  =  \g—  ny  [A^/n^J.  we  may  now 

i  =  l 

write  the  limiting  attacks  of  node  S^"  as 


and 


[Ag/ng  ]  ,  i  <  ig 

,  i — 1,  ...,  ns,  g — 1,  ...,  G. 

[Aj/n,  J,  I  >  ig 
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p 

With  ij  =  A,  —  the  limiting  attacks  of  node  S  are 

aR=a 

and 

f  fAf/nsli  1  <  It 

—  Ns  +  i  =  1  ’  ,=  ^t  •••>  nJt  9~  1 1  •••«  G. 

l  LAj/njJ,  «  >  Is 

3.2  Modification  of  the  Algorithm  for  Problem  97 

In  the  algorithm  of  section  2.2.2,  each  iteration  of  the  allocation  step  (step  Jb)  caused  only  the 
attack  at  the  “critical”  target  «’  to  be  changed.  In  the  modified  algorithm  below,  a  “critical”  group  is 
identified  with  an  allocation  to  all  targets  within  that  group.  This  faster  consumption  of  resources 
(attackers)  reduces  the  number  of  time-consuming  probability  calculations  that  needs  to  be  made.  We 
use  the  same  group  notation  as  introduced  above. 

Algorithm  for  Problem  with  Groups  (V>  convex): 

Step  0  (Initialization): 

a»  =  *=1»  T;  Aj  =  A; 

ik  =  1;  go  to  step  k. 

Step  k  (Allocation): 

If  At  =  0, 

stop  (a  is  optimal  for  problem  "?2); 

else 

x*  =  TAt/Tl; 
g*  = 

minj*  |  ®f[^(*Nf+i+x»)-^(«N,+i)]  =  +x* )  ~  ^(^^1  )]}•  9  =  1,  •••>  g|; 
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if  ntf*xt  <  Ak  then 

a.,  .  :=  a  .  4-  At, 

Nj  +i  Nj  +•  * 

At  +  i  :  =  Ak  —  n3mxk; 
k  k+  1;  go  to  step  k. 


I  —  If  ...»  n g  , 


a. |  .  :=  a  .  .  +  xt. 

Nj  +1  N  g  +  «  * 

a.  m  :  =  a. .  «  ,  4-  xt  1. 

Ng  +  1  Nj  +1  * 

stop  (a  is  optimal  for  ‘Jj). 


i=l,  Aj.  -  n,*(xt  -1); 

»  =  At  -  n,*(xt-l)4-l,  ...»  n,*; 


3.3  Modification  of  the  Algorithm  for  Problem 

When  Problem  *iP3  has  groups  of  equally  valued  targets,  we  may  modify  the  algorithm  of 
section  2.3.2  to  achieve  some  computational  efficiencies.  As  in  the  preceding  section,  allocation  of 
resources  at  each  step  now  applies  to  all  targets  within  a  group  rather  than  to  a  single  target.  Of 
course,  a  group  may  consist  of  only  one  target. 


Algorithm  for  Problem  ^3  with  Groups: 
Step  1  (Initial  Phase): 


a,  =  0,  i=l,  ...,  T; 


determine  Gf  =  jnaxc|j  |  £  n3\\  +  ln(vgl  j In  7J  <  a|; 

aNj+;  :=  U  +  T-l’  J=l,  ....  n(;  i  =  l,  ...,  G'-l; 

G'-i 

A  =  A  -  £  n<  aN  ; 

»  =  1  1 

aNG,  +  ,  =  ,=  1'  -1  min{A',  nQ,}; 

A1  :=  Ar  —  min{A,)  n  ,}; 

G 

if  A '  —  0,  stop  (a  is  optimal  for  <?3);  else  go  to  step  2. 
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Step  2  (Intermediate  Phase): 


if  A*  =  0,  stop  (a  is  optimal  for  93);  else  go  to  step  3. 


Step  3  (Final  Phase): 

Let  i,  iG,  be  a  reordering  of  1 . G*  such  that 

/  aN  + 1  aNi+i 

if  1  <  ]  <  k  <  G  ,  then  y  >  vky  ; 

determine  o'  =  max  <o  I  N,  <  AH; 

«=i . G'l  1  ~  J’ 

aN,?+h  =  aNtJ  +  h  +  A=l,...,n,s,  g=  1,  g’ 

A':=  A'-Ni#,; 

aN,s*+h  :=  *"ig*+h  +  h  *  = 
stop  (a  is  optimal  for  93). 


-  l; 


4.  COMPUTATIONAL  EXPERIENCE 

All  algorithms,  both  initial  and  modified,  have  been  implemented  in  FORTRAN  with  no 
special  compiler  options  invoked.  All  computations  have  been  performed  on  the  VAX  11/780  of  the 
School  of  Engineering  and  Applied  Science,  The  George  Washington  University. 


4. 1  Problem 

The  algorithm  for  problem  is  the  most  computationally  intense  of  the  three  in  that  it 
involves  two  levels  of  optimization  (both  attacker  and  defender),  each  of  which  involves  a  branch-and- 
bound  scheme.  Accordingly,  only  problems  of  modest  size  have  reasonable  run  times.  Table  1, 
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however,  shows  that  we  can  consider  solving  much  larger  problems  in  a  timely  manner  when  there  are 
groups  of  equally  valued  targets.  The  CPU  times  shown  in  the  “singly”  column  are  for  the  basic 
algorithm  in  which  targets  are  considered  individually.  The  times  shown  in  the  “Groups”  column  are 
for  the  modified  algorithm.  The  problem  data  are:  (T  =  )  20  targets,  (G  =  )  3  groups  (target  values: 
2@110,  12@60,  6@10),  attack  and  defense  reliabilities  a  —  6  =  0.9,  and  equal  numbers  of  attackers 
and  defenders  (A  =  D). 

Table  1.  Run  Times  for  Problem 


A(=D) 

CPU  Time  (secs) 

Reduction  Factor 

(Singly  -f  Groups) 

Singly 

Groups 

60 

844.03 

2.53 

333 

70 

1178.67 

3.30 

357 

80 

1637.69 

4.10 

399 

90 

1425.23 

4.75 

300 

100 

2091.03 

8.15 

256 

Total  Time 

7176.65 

22.83 

314 

The  significant  reduction  in  CPU  time  is  due  to  the  fact  that  the  attacker’s  branch-and-bound 
scheme  operates  in  G-space  in  the  modified  algorithm  rather  than  in  T -space,  as  in  the  basic 
algorithm. 

4.2  Problem 

For  the  fully  coordinated  defense  with  no  impact  point  prediction,  we  can  solve  larger  problems 
since  only  the  attacker  has  a  real  decision  on  allocating  resources.  In  Table  2  the  problem  data  are: 

(T  =  )  100  targets,  (G  =  )  3  groups  (target  values:  (10@22,  60@12,  30@2),  attack  and  defense 
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reliabilities  a  =  S  =  .99,  and  equal  numbers  of  attackers  and  defenders  (A  =  D). 


Table  2.  Run  Times  for  Problem 


A(  =  D) 

CPU  Time  (secs) 

Reduction  Factor 

(Singly  -r  Groups) 

Singly 

Groups 

200 

35.42 

.15 

236 

300 

112.94 

.28 

403 

400 

262.67 

.48 

547 

500 

506.78 

.82 

618 

600 

866.44 

1.30 

666 

Total  Time 

1784.25 

3.03 

588 

Here  the  reduction  in  run  time  can  be  attributed  to  the  much  smaller  number  of  survival 
probability  calculations  that  need  to  be  made.  This  is  a  direct  result  of  the  Quasiuniform  Attack 
Theorem. 

4.3  Problem  93 

The  exponential  structure  of  the  objective  function  in  the  model  formulation  for  Problem  ^P3 
permits  solution  of  problems  much  larger  than  either  of  the  others.  Notwithstanding  the  small  run 
times  of  the  basic  algorithm,  there  is  a  substantial  reduction  in  run  time  when  advantage  is  taken  of 
equally  valued  targets.  The  data  for  Table  3  are:  fixed  attack  size  of  (A  =  )  6000,  (T=)  100  targets, 
(G  =  )  3  groups  (target  values:  10@22,  60012,  3002),  attack  and  defense  reliabilities  a  =  .40, 

S  =  .80,  and  coordination  factor  (r=)  3. 
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Table  3.  Run  Times  for  Problem  ^P3 


D 

CPU  Time  (secs) 

Reduction  Factor 

(Singly  4-  Groups) 

Singly 

Groups 

14,000 

.31 

.02 

15 

15,000 

.27 

.02 

13 

16,000 

.31 

.02 

15 

17,000 

.29 

.02 

14 

18,000 

.29 

.02 

14 

Total  Time 

1.47 

.10 

14 

5.  CONCLUDING  REMARKS 

In  this  paper  we  have  emphasized  the  min-max  modeling  and,  especially,  the  solution 
efficiency  of  the  attacker/defender  problem.  We  have  done  so  for  three  different  cases  that  model 
alternate  degrees  of  defender  capability  with  respect  to  impact  point,  prediciton  and  coordination.  By 
tailoring  the  algorithms  to  the  situation  in  which  only  a  few  different  target  values  are  present,  we  are 
able  to  solve  problems  of  realistic  size  at  reasonable  cost.  Multiple  computer  runs  are  thus  feasible,  and 
so  it  becomes  practical  to  examine  sensitivity  and  trade-off  issues  by  solving  many  problem  variations. 

One  such  issue  we  have  already  begun  to  examine  is  the  sensitivity  of  results  to  the  attacker's 
perception  of  the  defender’s  capabilities  [see  O'Meara  and  Soland  (1988d)].  For  example,  an  attack 
that  is  optimal  against  a  defender  possessing  IPP  is  not,  in  general,  optimal  against  the  same  defender 
lacking  IPP.  With  efficient,  algorit  hms  for  the  three  cases  of  defensive  capability  examined  above,  it  is 
quite  easy  to  determine  the  penalties  paid  by  the  attacker  for  misjudging  the  capability  of  the  defense. 
The  attacker  optimization  models  developed,  and  the  defender  optimization  model  for  the  case  of  IPP. 
can  also  be  used  to  examine  the  sensitivity  of  results  to  other  parameters.  For  example,  the  attacker's 
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knowledge  of  a,  D,  6  and/or  r  may  be  imprecise.  Similarly,  the  defender’s  knowledge  of  a,  6  and/or  r 
may  be  imprecise.  Examination  of  the  effects  of  such  imprecision  may  be  carried  out  in  a 
straightforward  manner.  It  is  also  possible  to  compare,  as  a  function  of  the  true  state  of  affairs  and  the 
respective  states  perceived  by  attacker  and  defender,  the  actual  and  anticipated  expected  total  survival 
values. 
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